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ABSTRACT 


Lowland river degradation of agricultural areas is a problem not only in Poland but also 
in Europe and worldwide. The objectives of the present study were to determine the 
biodiversity of the gastropod communities in lowland rivers under the impact of agriculture, 
the relationship between gastropod communities and environmental factors, and to esti- 
mate the usefulness of the gastropod communities as an indicator of water quality. Four 
rivers that flow through an agricultural area were investigated. A high value of chlorophyll 
a, a low value of total dissolved oxygen and a high concentration of phosphates or total 
phosphorus and organic matters are present in the rivers. Gastropod communities in the 
rivers of the Ciechanowska Upland are influenced by bottom sediments, the physical and 
chemical parameters of water, and the abundance of macrophytes, which are typical of 
eutrophic water. In the rivers of the Ciechanowska Upland, gastropods can be biological 
indicators of dissolved oxygen in water — Viviparus viviparus (Linnaeus, 1758) and Lym- 
naea peregra (O. F. Miller, 1774) can be indicators of anthropopressure. The values of 
diversity indices varied from site to site, but did not decrease along the river. 

Key words: gastropod communities, lowland rivers, biodiversity indices, agricultural ar- 
eas, biological indicators, bottom sediments, Viviparus viviparus (Linnaeus, 1758), Lym- 


naea peregra (O. F. Miller, 1774). 


INTRODUCTION 


The cleanness of rivers reflects all human 
activities within a catchment (Hildrew, 1996). 
In the rivers of agricultural areas, high concen- 
trations of phosphates or nitrates accelerate the 
development of phytoplankton, diminishing 
macrophytes and increasing deficits in dis- 
solved oxygen (Thiébaut & Muller, 1999). It has 
been estimated that in inland waters, 70% of 
the total phosphorus comes from surface flow 
from agricultural areas, 24% comes from mu- 
nicipal sewage, and 19% from detergents (Parr 
& Mason, 2003). Most rivers in agricultural ar- 
eas have high concentration of total phospho- 
rus, phosphates, nitrates, suspended solids, 
and pesticides. This includes such rivers as the 
tydynia River, the Pełta River, the Sona River, 
and the Wkra River, which flow throughout the 
Ciechanowska Upland in central Poland. Low- 
land river degradation of agricultural areas is a 
problem not only in Poland but also in Europe 


(Solimini et al., 2000; Parr & Masson, 2003) 
and worldwide (Collier et al., 1998; Harding et 
al., 1999; Shieh et al., 1999). 

Biodiversity assessments are frequently 
used to monitor the level of pollution (Sandin 
& Jochnson, 2000; Beaven et al., 2001). Di- 
versity indices such as the Shannon-Wiener 
index and the Evenness index are based on 
the following assumptions: (1) there are fewer 
species in aquatic ecosystems under human 
impact compared to unpolluted water; (2) the 
diversity decreases with environmental deg- 
radation; and (3) values of the Shannon- 
Wiener index decrease with increasing 
pollution (Barton & Metcalfe-Smith, 1992). The 
objective of the present study were to deter- 
mine the biodiversity of the gastropod com- 
munities in lowland rivers under the impact of 
agriculture, the relationship between gastro- 
pod communities and environmental factors, 
and to estimate the usefulness of gastropod 
communities as indicators of water quality. 
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MATERIALS AND METHODS 
Study Area 


The study area comprised 2,570 km? (the 
Ciechanowska Upland). The study was carried 
out in the years 1997—2000. Four rivers of the 
Ciechanowska Upland: the tydynia River, the 
Pełta River, the Sona River and the Wkra River, 
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which flow through this agricultural area, were 
investgated. Samples were gathered from the 
overall length of the rivers, but in the case of 
the Wkra River, only the middle and lower 
course were sampled (103 km) (Fig. 1). Ten 
sampling sites were chosen for the Lydynia, 
Sona, Wkra, and eight for the Pelta River. 

The characteristic morphometric features of 
the rivers are shown in Table 1. 
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FIG. 1. Location of the study area (the Ciechanowska Upland). 
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TABLE 1. The morphometry of the rivers and the substratum. 


Catchment 


River length River width Surface velocity 


River area (inkm*) (in km) (inm) (ranges) (m/s) Substratum 
tydynia River 697.9 TES 1.5-12.0 0.17-0.91 sandy-stony, sandy, 
stony, stony-sandy 
Petta River 368.7 S0 0.7—4.0 0.03-0.3 sandy-muddy, sandy- 
stony, muddy-clay 
Sona River 536.5 67.7 0.5-15.0 0.10-0.98 stony, stony-sandy 
Wkra River 5322 255% 10.6-30.5 0.18-0.91 muddy-stony, sandy- 
muddy, sandy-stony 
Methods ried out by means of a Siemens D 5,000 pow- 


Gastropods were sampled by placing a 
quadrat frame (25 x 25 cm) on the substratum 
of the river. The frame was placed 16 times at 
each of the sample stations, which constituted 
one sample. At each sample station, samples 
were also taken of bottom sediments up to 5 
cm, macrophytes and water surfaces. Only liv- 
ing specimens of gastropds were collected. 
Bottom sediments were taken by a core-type 
sampler. The collected material was brought 
back to the laboratory in plastic bags. The 
samples were than filtered using a 0.5 mm 
mesh sieve. The samples of gastropods were 
preserved in 75% ethanol. The species of most 
gastropods were identified according to Gléer 
& Meier-Brook (1998), whereas the Lymnaei- 
dae were identified according to Jackiewicz 
(1998). The density of gastropods was esti- 
mated as the number of individuals per square 
metre. Immediately prior to gastropod sam- 
pling, water samples were collected from each 
sampling site. The analyses of the physical 
and chemical parameters of water, for ex- 
ample, temperature, pH, conductivity, dis- 
solved oxygen, BODs, sulphates, chlorides, 
total dissolved solids, alkalinity, calcium, mag- 
nesium, ammonia nitrogen, nitrite nitrogen, 
nitrate nitrogen, phosphates, total phospho- 
rus, and chlorphyll a, were carried out accord- 
ing to Hermanowicz et al. (1976). Macrophyte 
species were recorded on the same visit as 
gastropod sampling. If macrophyte species 
could not be identified to species in the field, 
they were taken to the laboratory, dried be- 
tween sheets of filter paper and after drying, 
mounted as ordinary herbarium species. Mac- 
rophytes were identified to species according 
to Szafer et al. (1986). The mineralogical 
analyses of the bottom sediments were car- 


der diffractometer. The organic matter content 
in the bottom sediments was determined ac- 
cording to Tiurin methods (Lityński et al., 
1976). All samples were collected from the 
head to the mouth of rivers three times at each 
sample site. In total, 114 sites were sampled 
within the survey area. The zoocenological 
study of the gastropod communities was car- 
ried out using the following indices: 
(1) Dominance (D%) 

D=n,/nx 100 

where n, = the number of individuals of 

species a, 

n = the total number of individuals in a 

sample. 

The value of the dominance index D was 
divided into five classes according to Gorny & 
Grim (1981): eudominants > 10.0% of the 
sample, dominants 5.1-10.0% of the sample, 
subdominants 2.1—5.05 of the sample, 
recedents 1.0-2.1% of the sample, sub- 
recedents < 1.0% of the sample. 

(2) Constancy (C %) 
C=N,/Nx 100 
where N, = number of samples that con- 
tain species a, 
N = total number of samples. 

The value of the constancy index C was di- 
vided into four classes: euconstants 75.1- 
100.0% of the sample, constants 50.1—75.0% 
of the sample, accessory species 25.1—50.0% 
of the sample, and accedents < 25.0% of the 
sample. 


Biodiversity indices (Cao et al., 1996; Mouillot & 
Lepretre, 1999): 
(1) The Shannon-Wiener index: 
m = -FP CRD 
where P; = N; / N — the proportion of indi- 
viduals belonging to species i. 
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(2) The Evenness index: 

J'=H’'/logs S$ 

where H’ = the value of the Shannon- 

Wiener index, 

S = the total number of species. 

Correlation between gastropod density, num- 

ber of species, organic matter content in the 
substratum, and the physical and chemical 
parameters of water were calculated by means 
of the non-parametric Spearman’s Rank Cor- 
relation Coefficient r,. Gastropod frequency in 
particular rivers in relation to the substratum 
was calculated by a chi-squared association 
test (v2) (Fovier et al., 1998). 


Principal Component Analysis 


The studied data sets contain measurement 
values of 18 various biological, physical and 
chemical parameters of water, collected at the 
sampling sites along the Pelta, tydynia, Sona 
and Wkra rivers. The water parameters stud- 
ied are listed in Table 2. The data are orga- 
nized in matrices X (35 x 18). Each row of 
matrix X represents one sampling site, de- 
scribed by 18 parameters. Because the mea- 
sured parameters significantly differ in their 
ranges, the data set was standardized accord- 
ing to the formula: 


where X,,5; denote the mean of the j-th col- 
umn and its standard deviation, respectively. 
Depending on the organization of the data 
sets, different methods of exploratory analy- 
sis can be used. One of the most popular tech- 
niques of exploratory analysis of multivariate 
data sets is principal component analysis 
(PCA) (Wold, 1987; Massart et al., 1997). This 
technique allows reduction of data dimension- 
ality, visualization and interpretation of the 
objects and variables relationships. 


RESULTS 


Physical and Chemical Parameters of Water 
and Macrophytes 


The physical and chemical parameters of 
water and macrophyte species are shown in 
Tables 2 and 3. 

A high value of chlorophyll a, a low value of 
total dissolved oxygen, and a high concentra- 
tion of phosphates or total phosphorus and 


TABLE 2. The physical and chemical parameters of water, organic matter content in the substratum, 


and their values. 


Parameters tydynia River 
Temperature (°C) 18.2-20.0 
pH 1./-8.2 
Conductivity (uS/cm) 490.0-610.0 
Dissolved oxygen (mg O2/dm°) 1.60-8.80 
BOD; (mg O2/dm°) 1.80-9.20 
Sulphates (mg SO4/dm°) 60.0-71.0 
Chlorides (mg Cl/dm°) 23.0-29.0 
Total dissolved solids (mg/dm*) 320.0-372.0 
Alkalinity (mg CaCO3/dm*) 200.0-245.0 
Calcium (mg Ca/dm*) 91.80-100.20 
Magnesium (mg Mg/dm°) 9.70-12.30 
Ammonia nitrogen (mg N/dm°) 0.30-3.74 
Nitrite nitrogen (mg N/dm”) 0.003—0.024 
Nitrate nitrogen (mg N/dm°) 0.072201 
Phosphates (mg PO4/dm°) 0.407210 
Total phosphorus (mg P/dm’) 0.48-1.01 
Chlorophyll a(ug/dm°) 1-80-2010 
Organic matter (%) 1.17—11.86 


Petta River Sona River Wkra River 
14.8—20.0 15.1-19.5 17.5-19.1 
7.6-7.9 7.6-8.4 8.0-8.5 
691.0—771.0 390.0-681.0 450.0-—620.0 
6.30-14.10 4.6—10.4 4.50-8.60 
1.30-7.80 1.20-5.60 1.20-8.0 
42.0-80.0 34.0-79.0 27.0-65.5 
20.0—51.0 15.0-23.0 19.0-22.0 
438.0-—613.0 280.0-488.0 295.0-369.0 
290.0-370.0 165.0-210.0 195.0—285.0 
90.70-104.70 64.90—104.0 77.0-101.1 
16.0—29.50 11.0-18.80 8.30—28.1 
0.86-1.98 0.06—1.16 0.27-0.92 
0.008-0.023 0.001—0.060 0.008—-0.095 
0.25-0.60 0.05-0.99 0.06-1.40 
0.34-3.71 0.30-0.60 0.68-1.31 
0.25-1.28 0.18-0.34 0.46-0.60 
3.90—55.1 7.40-81.50 3.8-75.6 

1.72-28.44 0.43-16.43 0.86-11.10 
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TABLE 3. Macrophyte occurrence in the rivers of the Ciechanowska Upland. 


Species 


Acorus calamus L. 

Alisma plantago-aquatica L. 

Berula erecta (Huds.) Coville 
Butomus umbellatus L. 
Ceratophyllum demersum L. s.s. 
Elodea canadensis Michx. 
Equisetum palustre L. 

Glyceria maxima (Hartm.) Holmb. 
Lemna minor L. 

Lemna trisulca L. 

Mentha agatica L. 

Myosotis palustris (L.) L. em. Rchb. 
Phragmites australis (Cav.) Trin. ex Steud. 
Nuphar lutea (L.) Sibth. & Sm. 
Potamogeton crispus L. 
Potamogeton natans L. 
Potamogeton pectinatus L. 
Potamogeton perfoliatus L. 

Rumex hydrolapathum Huds. 
Sagittaria sagittifolia L. 
Schoenoplectus lacustris (L.) Palla 
Sparganium emersum Rehmann 
Sparganium erectum L. em. Rchb. s.s. 
Spirodela polyrhiza (L.) Schleid. 
Typha latifolia L. 


2 of species 


organic matters were presentin the rivers. The 
main sources of pollution were public utilities, 
agriculture and the food industry. 

The organic matter in the bottom sediments 
ranged from 0.43% in the Sona River to 
28.44% in the Pełta River (Table 2). 


Zoocenological Study of the Gastropod Com- 
munities 


In the rivers of the Ciechanowska Upland, 
27 gastropod species were found (Table 4). In 
the Wkra River, the number of gastropod spe- 
cies was lower than in its tributaries, for ex- 
ample, the Lydynia River and the Sona River. 
The highest number of gastropod species were 
found in the tydynia River. Valvata piscinalis 
(O. F. Muller, 1774) was eudominant in the 
communities in the Lydynia River and in the 
Pełta River; Bithynia tentaculata (Linnaeus, 
1758) was eudominant in the Sona River, and 
Viviparus viviparus (Linnaeus, 1758) in the 
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Wkra River. Both Bithynia tentaculata and 
Lymnaea peregra (O. F. Müller, 1774) were 
eudominants and euconstants in communities 
in the Łydynia River, whereas Lymnaea 
peregra in the Wkra River. Seven gastropod 
species were subrecedents as well as 
accedents of the gastropod communities 
(Table 4). In the Lydynia River, the values of 
the Shannon-Wiener index and the Evenness 
index ranged from 0 to 3.28 and from 0 to 0.98, 
respectively. At sites 9 and 10, the values of 
H’ and J indices reached 0, because a drastic 
decrease in the number of gastropods spe- 
cies was observed there. The quantity ratio 
(the number of specimens of each species) 
was poorly balanced in the gastropod com- 
munity at site 5 (dominance was focused on 
the one species mainly, for example, Ancylus 
fluviatilis O. F. Muller, 1774), thus the value of 
the Evenness index reached 0.6 (Fig. 2a). In 
the Pelta River the values of the Shannon- 
Wiener and the Evenness indices ranged from 
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TABLE 4. The values of the dominance (D) and constancy (C) indices (%) of the gastropod commu- 


nities in the rivers. 


tydynia River Petta River Sona River Wkra River 

(n = 30) (n = 24) (n = 30) (n = 30) 

Species D C D C D C D C 

Viviparus contectus (Millet, 1813) 0.4 16.7 - - 3.4 20.0 - - 
Viparus viviparus (Linnaeus, 1758) - - Gye 42 12 Gyr 256 679 
Bithynia tentaculata (Linnaeus, 1758) 12.1 73.3 89 542 441 66.7 254 64.3 

Valvata cristata O. F. Müller, 1774 0.8 Of GO we C2 oo - - 

Valvata pulchella Studer, 1820 o7 J00 43 %0 0? 67 - - 
Valvata piscinalis (O. F. Müller, 1774) 224 40.0 27.7 29.22 1.0 20.0 46 28.6 
Valvata naticina Menke,1845 Or 293 - - - - 0.4 3:6 
Acroloxus lacustris (Linnaeus, 1758) 4.5 10.0 - - 02 33 ti 107 
Lymnaea stagnalis (Linnaeus, 1758) Ag “oeO 3:0) 583 57 63.3 3.7 -607 
Lymnaea palustris (O. F. Muller, 1774) 0.1 6.7 - - Or 23 03- 107 
Lymnaea corvus (Gmelin, 1791) 14 267 74 208 22 40.0 Ge 17.9 

Lymnaea truncatula (O. F. Muller, 1774) Gs 10 OT 42 04 10.0 - - 

Lymnaea auricularia (Linnaeus, 1758) - - G17 Zz 02 23 - - 
Lymnaea peregra (O. F. Müller, 1774) 12 ono 179 477 49 600 17.0 786 
Planorbis planorbis (Linnaeus,1758) E OSR S” Foe, 33.9. 22 214 

Planorbis carinatus O. F. Muller, 1774 0.1 33 - - - - - - 

Anisus spirorbis (Linnaeus, 1758) wa 167 70 458 04 16.7 - - 
Anisus vortex (Linnaeus, 1758) 62 393 3:0 438 36 36.7 3.4 50.0 
Bathyomphalus contortus (Linnaeus, 1758) 1.3 16.7 06 292 20 10.0 OFS. 36 
Gyraulus albus (O. F. Müller, 1774) 7 300" 02 12% 42 36.7 So) 290 

Gyraulus crista (Linnaeus, 1758) - - - - 0:2 33 - - 
Hippeutis complanatus (Linnaeus, 1758) - - - - - - On 5:10) 
Segmentina nitida (O. F. Muller, 1774) 0.1 33 - - G2 W00 02 107 
Planorbarius corneus (Linnaeus, 1758) ie Ao Se, 9000 96 7335 44 393 
Ancylus fluviatilis O. F. Miller, 1774 182 300 Of 42 - - AD 32 
Physa fontinalis (Linnaeus, 1758) Ge “der 42 730 56 500 4.7 57.1 

Aplexa hypnorum (Linnaeus, 1758) G2 133 12 22 .- - - - 

2 of species 23 19 22 18 
2 of specimens 2,228 2,794 1977 2,654 


0.0 to 2.80 and from 0.0 to 0.88, respectively 
(Fig. 2b). From site 5 to site 8, the values of H’ 
indices showed a decreasing and increasing 
tendency; at site 5 the value of H’ index 
reached 1.98, at site 6 the value of H’ reached 
2.71, at site 7 and 8 values reached 1.94 and 
2.32, respectively. The Evenness index value 
increased from site 5, where the structure of 
gastropod dominance started equalizing. In the 
Sona River the values of H’ and J indices 
ranged from 0.81 (site 1) to 2.93 (site 10), and 
0.43 (site 8) to 0.91 (site 4). The maximum 
value of H’ = 2.93 for the gastropods commu- 
nity occurred at site 10 (Fig. 2c). The values 
of H’ and J indices ranged from 1.4 to 2.83 or 
0.61 to 0.71, respectively, in the Wkra River 
(Fig. 2d). 


Analysis of the Spearman Rank Correlation 
Coefficent r, 


The Spearman Rank Correlation Coefficent 
r, (n = 38) showed a lack of statistically sig- 
nificant correlations between gastropod den- 
sity (r, = 0.36, P= 0.12), the number of species 
(r, = 0.06, P = 0.80) and the organic matter 
content in the bottom sediments. There were 
statistically significant positive correlations be- 
tween gastropod density and the number of 
species (r, = 0.65, P< 0.001) and magnesium 
concentration (r, = 0.44, P < 0.05) in the wa- 
ter, the number of species and dissolved oxy- 
gen concentration (r, = 0.46, P < 0.05) (Table 
5). The density of gastropods varied from 10 
to 135 individuals/m2 below 11.0 mg Mg/dms, 
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FIG. 2. The value of H’ and J for the gastropod communities in the rivers of the Ciechanowska 
Upland. (a): The value of H’ and J for the gastropod communities in the Lydynia River; (b): The value 
of H’ and J for the gastropod communities in the Petta River; (c): The value of H’ and J for the 
gastropod communities in the Sona River; (d): The value of H’ and J for the gastropod communities 
in the Wkra River. 


TABLE 5. The values of the Spearman Rank Correlation Coefficent rs. 


Value of the Spearman Rank 


Pair of factors Correlation Coefficent rs 
Specimen density and number of species 0:66 P= 0,001 
Specimen density and magnesium concentration 0.44 P<0.05 
Number of species and dissolved oxygen concentration 0.46 P<0.05 
Density of Viviparus viviparus and pH 0.51 PSA 
Density of Viviparus viviparus and BODs Cael we Os 
Density of Viviparus viviparus and chlorophyll a concentration 0.38 "P= 001 
Density of Lymnaea peregra and conductivity -0.54 P<0.01 
Density of Lymnaea peregra and calcium concentration -0.47 P<0.05 


Density of Lymnaea peregra and sulphates concentration -0.47 P<0.05 
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TABLE 6. Gastropod distribution depending on the parameters of the water. 


Parameter Range Density (individuals/m’) 
Density of Viviparus viviparus depending on pH, BODs and chlorophyll a 
pH < 8:0 0-36 
8.1-8.3 0-70 
> 8.4 17-96 
BOD; (mg O2/dm*) $38 0-38 
3.9-5.0 0-70 
Gr 21-74 
Chlorophyll a (ug/dm’) < 10.0 0-36 
10.1-15.0 3—39 
Salton 1—96 
Density of Lymnaea peregra depending on conductivity, calcium and sulphates 
Conductivity (uSxcm’) < 500.0 12-90 
501.0-—600.0 0-34 
> 601.0 0-18 
Calcium (mg Ca/dm*) < 85.0 1-90 
85.1—95.0 3-61 
= Qa: 0-18 
Sulphates (mg SO4/dm°) < 40.0 12-90 
40.1-60.0 1-61 
> 60.1 0-13 


from 10 to 223 individuals/m2 at 11.1-17.8 mg 
Mg/dm? and from 36 to 230 individuals/m2 at 
above 17.9 mg Mg/dm3. The number of gas- 
tropod species varied from 1 to 9 below 6.6 
mg O,/dms in the water, from 3 to 9 at 6.7— 
7.1 mg O.,/dm? and from 6 to 10 at above 7.2 
mg O./dms. There were statistically significant 
positive correlations between the density of 
Viviparus viviparus and pH, BOD,, chlorophyll 
a, and a negative correlation between den- 
sity of Lymnaea peregra and any of conduc- 
tivity, calcium and sulphate concentration 
(Tables 5, 6). 


Fequency Analysis of Gastropod Communities 


The result of the frequency analysis showed 
statistically significant positive (+), negative (-) 
or lack of associations (0) between certain gas- 
tropod species, for example, Bithynia tentacu- 
lata, Lymnaea stagnalis (Linnaeus, 1758), 
Lymnaea peregra or Anisus vortex (Linnaeus, 
1758) and the substratum in the Lydynia River 
(x2, = 64.46, P < 0.01), the Petta River (x2, = 
205.64, P< 0.01), the Sona River (72, = 62.03, 
P < 0.01), and the Wkra River (72, = 175.8, P< 
0.01). Bithynia tentaculata was positively as- 


TABLE 7. Gastropod frequency on the substratum in the rivers of the Ciechanowska Upland. 


Species tydynia River 


Bithynia tentaculata — sandy-stony (+) 
Lymnaea stagnalis 
Lymnaea peregra 
Anisus vortex 


sandy (-) 
sandy (+) 
sandy-stony (+) 
muddy (-) 


Peita River 


muddy-clay (-) 


muddy-clay (-) 
muddy-stony (-) 
muddy-stony (-) 


Sona River Wkra River 


sandy-muddy (-) 
sandy-stony (+) 
muddy (+) 
muddy (+) 
sandy-muddy (+) 


sandy-stony (0) 


sandy-stony (0) 
sandy-stony (0) 
sandy-stony (+) 
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sociated with a sandy-stony substratum in the 


tydynia River and in the Wkra River an avoid- 


ing muddy-clay or sandy-muddy substratum. 
In the Wkra River, Lymnaea stagnalis was posi- 
tively associated with a muddy substratum. 
Lymnaea peregra was positively associated 
with a sandy and muddy substratum, whereas 
Anisus vortex was positively associated with a 
sandy-stony (Lydynia River) and sandy-muddy 
substratum (Wkra River) (Table 7). 


Principal Component Analysis Results 
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To explore the studied monitoring data set and 
to examine the similarities of the samples, a 
chemometrics technique of the Principal Com- 
ponent Analysis was used. As this data set con- 
tains measurements performed within different 
magnitude ranges, the PCA model was con- 
structed for the centered and standardized data. 
Unfortunately, the reduction of data dimen- 
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FIG. 3a. Score and loading plots as a result of PCA for centered and standardized data X (35 x 18). 
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FIG. 3b. Score and loading plots as a result of PCA for centered and standardized data X (35 x 18). 


sionality was not too effective, because the PCA 
model with seven significant Principal Compo- 
nents describes 87.3% of data variance. The 
score plots and loading plots, which were ob- 
tained as a result of this analysis, are presented 
in Figs. 3a, D. 

Based on Figs. 3a, b it is possible to con- 
clude that PC1 reflects the difference between 
sample sites 8—10 (Sona River), 1-4, 6, 7 (Wkra 
River), and all remaining samples. Moreover, 
PC1 reveals the uniqueness of sample site 2 
(Petta River). Based on the loading plot PC1- 


PC2, it can be stated that sample sites 8—10 
(Sona River), 1—4, 6, 7 (Wkra River) differ from 
the remaining ones, mainly due to the high val- 
ues of variables 2 (number of species), 3 (pH), 
13 (nitrite nitrogen), 17 (chlorophyll a) and 14 
(nitrate nitrogen). Moreover, PC1 shows the 
uniqueness of sample site 2 (Pełta River) due 
to the highest value of parameter 7 (chlorides). 

PC2 is constructed mostly due to the differ- 
ence between sample site 5 (Sona River) and 
sample sites 2, 3 (Petta River) 8 and 6 (Lydynia 
River). Sample site 5 (Sona River) is charac- 
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terized by a high value of parameter 5 (dis- 
solved oxygen), whereas sample sites 2, 3 
(Petta River) 8 and 6 (Lydynia River) are char- 
acterized by the highest value of parameters 6 
(BOD,), 12 (ammonia nitrogen) and 15 (phos- 
phates), as well as the relatively lower values 
of the remaining parameters. | 

PC3 shows the difference between sample 
site 3 (Petta River) and sample sites 7 and 8 
(Lydynia River). Sample site 3 (Petta River) is 
characterized by a high value of parameters 1 
(gastropod density), 6 (BODs), 11 (magnesium), 
13 (nitrite nitrogen) and 14 (nitrate nitrogen), 
whereas sample sites 7 and 8 (Lydynia River) 
are characterized by the highest value of pa- 
rameter 17 (chlorophyll a). 

PC4 reveals the uniqueness of sample site 1 
(Pelta River), which is characterized by the low- 
est value of parameter 2 (number of species), 
whereas PC5 shows the uniqueness of sample 
sites 5 and 6 (Pelta River) characterized by the 
highest value of parameter 18 (density of Lym- 
naea peregra). 

Furthermore, PC6 is constructed mainly due 
to the difference between sample site 8 and 
sample site 1 (Pelta River), whereas PC7 is 
constructed due to the difference between 
sample sites 8 (Sona River) and 5 (tydynia 
River). Sample site 8 (Pelta River) is charac- 
terized by the highest value of parameter 13 
(nitrite nitrogen), whereas sample site 1 (Pelta 
River) is characterized by the highest value of 
parameter 1 and 9 (gastropod density and al- 
kalinity). Based on PC7 it is possible to con- 
clude that sample site 8 (Sona River) is 
characterized by the highest values of param- 
eter 17 (chlorophyll a), whereas sample site 5 
(Lydynia River) is characterized by the high 
values of parameters 14 (nitrate nitrogen), 16 
(total phosphorus), 18 (density of Lymnaea 
peregra) and 5 (dissolved oxygen). 

The PCA loading plots are easy to interpret 
and reveal a high positive correlation between 
the parameters, namely: 

— cluster 1: parameters 10 and 8 (calcium and 
sulphates), 

— cluster 2: parameters 4 and 9 (conductivity 
and alkalinity), 

— cluster 3: parameters 11 and 7 (magnesium 
and chlorides), 

— cluster 4: parameters 15 and 16 (phosphates 
and total phosphorus), and 

— cluster 5: parameters 6 and 12 (BOD; and 
ammonia nitrogen). 

Based on the loading plots it is also possible 
to observe the negative correlation between 
parameters 17 and 5 (chlorophyll a and dis- 
solved oxygen). 


DISCUSSION 


Gastropod Communities and their Relation with 
the Basement Complex and the Substratum 


Data of the surveys conducted by Harding 
et al. (1999) or Dietrich & Anderson (2000) 
showed that basement complex, substratum, 
and current velocity influenced the gastropod 
communities as well as the organic matter 
content in aquatic ecosystems (Vannote et al., 
1980; Delong & Brusven, 1998; Murphy & 
Giller, 2000). 

According to Savage & Pratt (1976), the 
percentage of organic matter content in the 
bottom sediments influences the macrophytes 
richness. In the rivers of the Ciechanowska 
Upland, the organic matter content in the bot- 
tom sediments was low and reached 0.43- 
28.44%. | did not observe any correlation 
between gastropod density and either the 
number of species or the organic matter con- 
tent of the water. This result is contrary to that 
of Savage & Gazey (1987), who obtained a 
statistically significant correlation between 
Bithynia tentaculata occurrence and the con- 
tent of organic matter. The data from the sur- 
vey of Dillon (2000), which involved both a 
substratum and a habitat size, showed an in- 
direct effect of calcium on the gastropod com- 
munities. In the New River and its tributaries 
investigated by Dillon, which passed through 
a region of gneiss and schist and then through 
limestone and dolomite, calcium did not affect 
the abundance of pulmonates directly. The 
gastropods were rare independent of the cal- 
cium content in the water. According to Dillon, 
the gastropods were also rare in hard water 
where the tributaries were small. Rosaro & 
Pietrangelo (1993) claimed that among 
macroinvertebrates, the gastropods are domi- 
nants in streams with a limestone bottom. The 
result of my survey showed the occurrence 
mainly of calcite (CaCO3) and dolomite 
(CaMg[CO3]>). These minerals could easily 
dissolve in water. The calcium content in the 
water varied from 64.90 to 104.70 mg/dms. 
Thus, such concentration of calcium could re- 
sult in a statistically significant negative cor- 
relation between Lymnaea peregra density 
and the calcium content in the water. Collier 
et al. (1998) came to the conclusion that the 
gastropod density is negatively correlated with 
a muddy and a sandy substratum. Strzelec 
(1993) also showed the negative influence of 
a muddy substratum on gastropod richness, 
where only Lymnaea stagnalis and Lymnaea 
peregra were present. The results | obtained 
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confirmed their surveys. In the tydynia River, 
gastropod richness was poor on a muddy sub- 
stratum. Besides, | found statistically signifi- 
cant negative association of Anisus vortex with 
a muddy substratum and Lymnaea stagnalis 
ona sandy substratum. But in the Wkra River, 
Lymnaea stagnalis and Lymnaea peregra 
were positively associated with a muddy sub- 
stratum. Bithynia tentaculata in the rivers of 
the Ciechanowska Upland were associated 
with a sandy-stony substratum but avoid a 
muddy-clay and sandy-muddy substratum. 
This result confirmed Dussart’s (1979) results 
that Bithynia tentaculata occurs on the sedi- 
ment surface and to a lesser degree is asso- 
ciated with macrophytes. The data of surveys 
done by Chertoprud & Udalov (1996) showed, 
that in the rivers 19-20 gastropod species 
occurred on macrophytes, 12—13 species on 
a substratum consisting of macrophytes and 
mud and detritus, and 7—10 species on a stony 
substratum. By contrast, my survey showed, 
that particular gastropod species were asso- 
ciated with a particular substratum. 
Stanczykowska (1960a) and Piechocki (1969) 
found gastropod diversity and abundance were 
dependent on the specificity of the study area. 
In the rivers of the southern part of Poland, for 
example, Anisus spirorbis (Linnaeus, 1768), 
Acroloxus lacustris (Linnaeus, 1758), and spe- 
cies of the genus Valvata were totally absent 
(Strzelec, 1993). These species were detected 
in the rivers of the Ciechanowska Upland. 
According to Hubendick (1964, 1972), 
Acroloxus lacustris is typical of reservoirs. By 
contrast, Acroloxus lacustris was subdominant 
in the Lydynia River, subrecedent in the Sona 
River and recedent in the community in the 
Wkra River. 


Gastropod Communities and their Relation to 
Physical and Chemical Parameters of Water 


According to Strzelec (1993), the number of 
gastropod species is negatively correlated with 
the sulphate and magnesium concentration in 
water if the sulphate concentration level 
reached above 210.0 mg SO,/dm3. | have 
found that magnesium concentration is posi- 
tively correlated with gastropod density in the 
rivers of the Ciechanowska Upland. Gastro- 
pod density was highest when magnesium 
concentration reached a level of 17.9 mg Mg/ 
dm’. The density of gastropods also de- 
creased when the magnesium concentration 
decreased. My result is consistent with the in- 
vestigation of Shieh et al. (1999). They found 


the highest gastropod density in rivers when 
the magnesium concentration ranged from 
14.55 to 77.27 mg Mg/dms. Dillon (2000) 
showed that the addition of magnesium to 
water improved the performance of freshwa- 
ter snails, provided that the concentration of 
magnesium did not exceed the concentration 
of calcium. Magnesium may compete for cal- 
cium uptake sites. In the rivers of the 
Ciechanowska Upland, the concentration of 
magnesium did not exceed the concentration 
of calcium. 

In the rivers of the Ciechanowska Upland, 
the sulphate concentration influences Lym- 
naea peregra density. The result of my inves- 
tigation confirmed the survey done by 
Clenaghan et al. (1998) that showed a posi- 
tive statistically significant correlation between 
gastropod density and magnesium concentra- 
tion as well as pH. They also found a positive 
correlation between gastropod density and 
conductivity, density, and calcium concentra- 
tion, whereas the results of my survey showed 
a negative correlation. 

Dissolved oxygen seems to play a role in 
gastropod density (Mouthon, 1996). Boycott 
(1936) and Hubendick (1972) found that 
Ancylus fluviatilis demanded water with a high 
level of dissolved oxygen. According to Collier 
et al. (1998), gastropod density is positively 
correlated with dissolved oxygen. In the riv- 
ers of the Ciechanowska Upland, there was a 
correlation between gastropod density and the 
concentration of dissolved oxygen. The num- 
ber of gastropod species varied from 6 to 10 
at above 7.2 mg O2/dm3 and from 1 to 10 be- 
low 7.2 mg O,/dm3. The results of my survey 
are consistent with the opinions of these au- 
thors. 

Dillon (2000) found that the diversity of fresh- 
water gastropods was at a maximum in wa- 
ters of about 5-40 mg Ca/dm3. Dussart (1976) 
claimed that the occurrence of Bithynia 
tentaculata is limited to hard- and medium-hard 
water. The calcium content in the rivers of the 
Ciechanowska Upland varied from 64.90 to 
104.70 mg Ca/dm3 and influenced Lymnaea 
peregra density only. If the calcium concen- 
tration increased up to 95.1 mg Ca/dms, the 
density of Lymnaea peregra decreased. If the 
calcium concentration reached 85.0 mg Ca/ 
dms, the density of Lymnaea peregra in- 
creased. Thus, the results of my survey con- 
firmed the fact that the number of specimens 
per m2 reaches higher levels in hard water 
according to Dussart’s classification. Boycott 
(1936) showed that Lymnaea peregra was 
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most common in water with a calcium concen- 
tration from 1 mg to 8 mg Ca/dm3. Briers 
(2003) found that variation in the environmen- 
tal calcium requirement is an important factor 
in determining the geographical ranges of 
freshwater gastropods. According to the envi- 
ronmental requirements, Lymnaea peregra is 
classified as a non-calciphilous species (Boy- 
cott, 1936; Briers, 2003). According to 
McMahon (1983), more than 45% of freshwa- 
ter basommatophoran snails occur only in rela- 
tively hard waters, that is, above 25 mg Ca/ 
dm°. For example, the optimal value of cal- 
cium concentration in water for Lymnaea 
palustris (O. F. Muller, 1774) amounts to 30.0- 
280.0 mg Ca/dm?, whereas for Ferrissia 
rivularis (Say, 1817) it is 4.6-67.6 mg Ca/dm%. 
Growth rate, survivorship, or fecundity rates 
all decline with a reduction in calcium concen- 
tration below optimal values. In the rivers of 
the Ciechanowska Upland, calcium concen- 
tration from 40.0 to 85.0 mg Ca/dm* may play 
a role in gastropod density. Above this con- 
centration, gastropod density decreased. Ac- 
cording to Young (1975), if the level of calcium 
in water is higher than 5.0 mg Ca/dm‘°, other 
factors determine the abundance and density 
of gastropods. The negative correlation be- 
tween Lymnaea peregra density and a high 
calcium content in water may explain the sta- 
tistically significant association of this species 
with a muddy substratum. A muddy substra- 
tum reduces calcium concentration in water 
by colloid absorption. Calcium does not accu- 
mulate in water until a muddy substratum is 
saturated with calcium (Dussart, 1976). Thus, 
a muddy substratum “protects” Lymnaea 
peregra against an excessive calcium concen- 
tration in water. Shieh et al. (1999) obtained a 
negative correlation between alkalinity, con- 
ductivity and gastropod richness, as well as a 
positive one between alkalinity and gastropod 
density. In the rivers of the Ciechanowska 
Upland, conductivity negatively influences 
Lymnaea peregra density only. If the conduc- 
tivity increases above 601.0 uS x cm”, the 
density of Lymnaea peregra decreases. | have 
found a statistically significant correlation be- 
tween conductivity and gastropod density, but 
not between conductivity and the number of 
gastropod species. 

Jokinen (1985) claimed that the number of 
the pulmonate species is correlated with al- 
kalinity and the catchment area. This result of 
the survey is easy to understand, because if 
the level of alkalinity is high, the rate of the 
allochtonous organic matter exchange is high 


compared to the lower value of alkalinity in 
smaller rivers. Thus, rivers with a higher value 
of alkalinity can feed more invertebrate spe- 
cies. Clarke & Scruton (1997) claimed that in 
rivers the absence of gastropods is connected 
with the pH values of water. This conclusion 
confirmed an investigation by Herrmann et al. 
(1993), who noticed a rapid decrease in the 
number of gastropod species if the pH de- 
creased to 5.5. In the rivers of the Cie- 
chanowska Upland, the density of Viviparus 
viviparus decreased if the pH value decreased. 
The Ciechanowska Upland is a typical agri- 
cultural area. The physical and chemical analy- 
sis of the water showed high concentration of 
nitrate nitrogen, total phosphorus, and phos- 
phates. Soil fertilization and discharge of mu- 
nicipal and industrial sewage affect the level 
of chlorophyll a in the rivers of the Ciecha- 
nowska Upland. | have obtained a statistically 
significant correlation between Viviparus 
viviparus density and the chlorophyll a level. 
There is no correlation between the density of 
Viviparus viviparus and either nitrate nitrogen 
or phosphates. However, the level of chloro- 
phyll a in water is dependent on the concen- 
tration of these nutrients. Thus, Viviparus 
viviparus may play the role of a bioindicator of 
eutrophication in running waters. 


Gastropod Communities and their Relation to 
Macrophytes 


Some gastropod species are associated with 
bottom sediments as well as with macro- 
phytes, for example, Bithynia tentaculata 
(Kornijow, 1989). According to Soszka (1975) 
Bithynia tentaculata and Valvata piscinalis 
occur on bottom sediments mainly, and to a 
lesser degree on macrophytes. Anisus vortex 
is associated mainly with bottom sediments 
(Kornijow et al., 1990). 

In the rivers of the Ciechanowska Upland, 
gastropods were frequently observed on mac- 
rophytes. In the Lydynia River, at one sam- 
pling site, Ancylus fluviatilis occurred 
exclusively on Nuphar lutea (L.) Sibth. & Sm. 
Planorbis carinatus O. F. Muller, 1774, a spe- 
cies typical of lakes, was present on Pota- 
mogeton crispus L. Stanczykowska (1960b) 
and Dvorak & Best (1982) suggested that 
Bithynia tentaculata, Lymnaea stagnalis, Lym- 
naea peregra, and Anisus vortex or Gyraulus 
albus (O. F. Muller, 1774) were not typical 
phytophilous species. | found that these spe- 
cies occurred on macrophyte surfaces as well 
as on bottom sediments. This result confirmed 
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the survey of Kornijow & Gulati (1992). Brown 
(1997) claimed that gastropods prefer wider 
the leaf blades of macrophytes rather than 
narrow leaf blades because they provide more 
periphyton. According to Brénmark (1989), 
Ancylus fluviatilis occurs on the lateral sur- 
faces of stones because they are richer in dia- 
toms. Valvata piscinalis feeds mainly on 
detritus, whereas green algae and diatoms 
constitute only a small part of their diet 
(Kornijow, 1996). McMahon (1983) showed 
that bottom sediments play a role for freshwa- 
ter gastropods because minerals comprise a 
large component of the periphyton on which 
many gastropods graze. Calcium in the shells 
derives both from ingested food and bottom 
sediments. Thomas & Kowalczyk (1997) 
claimed that the relationship between macro- 
phytes, periphyton, and gastropods is based 
on mutualism, whereas Jones et al. (2000) 
obtained the opposite result. According to 
Dillon (2003), the effects of calcium, alkalinity, 
conductivity or pH on gastropod communities 
may be indirect. Thus, higher values of these 
parameters may influence the production of 
food on which gastropods might feed. 
Kołodziejczyk (1984) found that Bithynia 
tentaculata and Anisus spirorbis prefer detri- 
tus. Thus, it is likely, that the statistically sig- 
nificant association of those species with the 
bottom sediments in the rivers of the 
Ciechanowska Upland may be explained by 
alimentation relationships. 


Gastropods, Macrophytes, and their Bio- 
indicative Values 


Many authors consider gastropods as a good 
indicator of sewage, heavy metals or radionu- 
clide concentration in rivers (Dregolskaya, 
1993; Abdallah et al., 1999; Flessas et al., 
2000). Shieh et al. (1999) claimed that gas- 
tropods could be indicators of water hardness. 
Gastropods are more sensitive to pH changes 
than fish (Okland, 1983). Calcium is not eas- 
ily accessible to gastropods at a low pH value 
(Clarke & Scruton, 1997). Thus, gastropods 
are good indicators of the acidification of wa- 
ter environments. The results of my survey 
show that certain physical and chemical pa- 
rameters influenced the optimal gastropod 
population development. According to 
Mouthon (1996), Valvata piscinalis and 
Viviparus viviparus showed an association with 
aquatic environments that are abundant in 
phosphates and nitrates. In the rivers of the 
Ciechanowska Upland, the correlation be- 


tween Viviparus viviparus density and the 
BOD; level was statistically significant, so that 
a moderate load of organic matter appeared 
to promote its occurrence. Valvata piscinalis 
and Viviparus viviparus were numerous in the 
middle and lower course of rivers, where the 
concentration of phosphates and total phos- 
phorus was high. In the rivers of the 
Ciechanowska Upland, gastropods could also 
be indicators of the concentration of dissolved 
oxygen. The number of gastropod species in- 
creased with an increasing level of dissolved 
oxygen in water. Viviparus viviparus and Lym- 
naea peregra could be indicators of anthropo- 
pressure. Viviparus viviparus density increased 
when pH, BOD; and the level of nutrients also 
increased. 

According to Thiébaut & Muller (1999), the 
occurrence of Elodea canadensis Michx. in the 
middle course indicates a progression of 
eutrophication, whereas Berula erecta (Huds.) 
Coville or Lemna trisulca L. indicates olig- 
otrophic water. In the case of my survey, Elo- 
dea canadensis occurred in the middle course 
of the rivers, but Berula erecta occurred in the 
upper course of the Petta River only. 


Diversity Indices 


According to Trojan (2000), gastropod com- 
munity diversity depends on the size of the 
study area, basement complex, and 
anthropopressure. Clausen & Biggs (1997) 
showed that the values of the Shannon-Wiener 
index are high for invertebrate communities 
and among these for gastropod communities 
that occurred in smaller rivers compared to 
communites in larger ones. | have obtained a 
similar result for the gastropod communities 
which occur in the Lydynia River and the Sona 
River, which are tributaries of the Wkra River. 
Chedwick et al. (1986) found that the Shan- 
non-Wiener index values decreased in rivers 
under the influence of the metal and mining 
industries. By contrast, Garcia-Criado et al. 
(1999) claimed that the Shannon-Wiener in- 
dex does not reflect the environmental impact 
of coal mines on community structure in riv- 
ers. Perhaps, if they had carried out their sur- 
vey at genus level instead of the family level, 
the environmental influences would have been 
more distinctive. In my survey, the effects on 
physical and chemical parameters differed 
from effects caused by the mining industry. The 
values of H’ and J indices were high for gas- 
tropod communities in the headwater of the 
tydynia River only. In the lower course of the 
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tydynia River, the values decrease to 0. A simi- 
lar result was obtained by Shieh et al. (1999) 
in rivers under the influence of industry and 
agriculture. Strzelec et al. (1999) also claimed 
that along rivers of the Sląska Upland (south- 
ern Poland) gastropod communities decrease. 
| have not observed a similar tendency in the 
Sona, Pełta, and Wkra rivers. In the rivers of 
the Ciechanowska Upland, diversity indices 
varied from site to site, but did not decrease 
along the river pollution gradient. The similar 
tendency was observed by Cao et al. (1996) 
in the River Trent system: diversity indices 
varied widely between sites, but did not show 
a consistent decrease with increasing pollu- 
tion levels; low values were recorded at both 
clear and polluted sites. 


CONCLUSIONS 


Gastropod communities in the rivers of the 
Ciechanowska Upland are influenced by bot- 
tom sediments, the physical and chemical 
parameters of water, and macrophyte abun- 
dance, which are mainly typical of eutrophic 
water. Valvata naticina Menke, 1845, threat- 
ened with extinction in Poland, is present in 
the Łydynia River and in the Wkra River. The 
result of the survey showed that the range of 
certain physical and chemical parameters in- 
fluenced the development of optimal gastro- 
pod population and the threshold values at 
which inhibition of their development occurred. 
In the rivers of the Ciechanowska Upland, 
gastropods can be biological indicators of the 
dissolved oxygen level in water. When the level 
of dissolved oxygen increased, the number of 
gastropods species also increased. Viviparus 
viviparus and Lymnaea peregra can be indi- 
cators of anthropopressure. Viviparus 
viviparus density increased with increases in 
pH, BODs, and the concentration of chloro- 
phyll a. Lymnaea peregra density decreased 
if calcium, sulphates, and conductivity in- 
creased. The values of diversity indices var- 
ied from site to site, but did not decrease along 
the river pollution gradient, so the result | have 
obtained is not easy to interpret. Perhaps the 
diversity indices do not reflect the 
anthropopressure degree in the rivers of the 
Ciechanowska Upland. Diversity indices may 
be more useful for aquatic ecosystems that 
are more polluted than the rivers of the 
Ciechanowska Upland. The PCA loading plots 
mainly showed a high correlation between the 
physical and chemical parameters. 
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